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Abstract 

In an influential paper of 2002, Babaev, Faddeev and Niemi conjectured that two- 
component Ginzburg-Landau (TCGL) theory in three dimensions should support knot 
solitons, where the projective equivalence class of the pair of complex condensate fields 
[^1,^2] : K 3 CP 1 has non-zero Hopf degree. The conjecture was motivated by a 
certain truncation of the TCGL model which reduced it to the Faddeev-Skyrme model, 
long known to support knot solitons. Physically, the truncation amounts to ignoring the 
coupling between [^1,^2] and the supercurrent of the condensates. The current paper 
presents a direct test of the validity of this truncation by numerically tracking the knot 
solitons as the supercurrent coupling is turned back on. It is found that the knot solitons 
shrink and disappear as the true TCGL model is reached. This undermines the reasoning 
underlying the conjecture and, when combined with other negative numerical studies, 
suggests the conjecture, in its original form, is very unlikely to be true. 

PACS: 11.27+d, ll.lOLm, 05.45Yv 



1 Introduction 

Two component Ginzburg-Landau (TCGL) theory is a phenomenological field theory which 
is enjoying increasing prominence in several areas of condensed matter theory. It consists of 
two complex scalar fields ipi, ip 2 ' K 3 —> C minimally coupled to a £7(1) gauge field A 6 ^ 1 (]R 3 ) 
and (perhaps) each other via an energy density of the form 

€ = \\d A H 2 + \\d A H 2 + \\dA\ 2 + U(^M (1) 

where d^ a = di/j a — iA\\) a and U is a potential function whose details depend strongly on the 
precise physical context. In an influential and much-cited paper of 2002 [2], Babaev, Faddeev 
and Niemi conjectured that models of this type should possess "knot solitons" in which the 
projective equivalence class of the pair of complex fields ^1, ^2, 

[^2] :M 3 ^CP 1 * S\ (2) 
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has non-zero Hopf degree. This conjecture is based on the observation that, when expressed 
in terms of the gauge-invariant fields p — ^\ipi\ 2 + |^2p and 

2 

C = ^ ^2(^ a d A ^ a - ^ a d A ^a) (3) 
a=l 

the TCGL energy density is 

£ = \p 2 \M 2 + \\dc + \vM 2 + \\M 2 + \p 2 \c\ 2 + u( P , <p), (4) 

where uu is the area form of the usual metric on the unit two-sphere. In superconductivity 
contexts, p 2 C is the total supercurrent associated with condensates ^1^2 and so, in a slight 
abuse of terminology, we shall refer to C itself as the supercurrent. If we impose that p = 
Po > 0, a constant, and C = 0, £ coincides with the energy density of the Faddeev-Skyrme 
(FS) model, which is known to possess knot solitons of every Hopf degree Q. These are 
smooth, spatially localized global energy minimizers in their homotopy class. Topologically, 
maps Lp : M 3 U {00} — >• S 2 are classified by the framed cobordism class of p)~ x {p) where p e S 2 
is any regular value of cp. In general, p>~ l {p) is a union of closed curves in M 3 , that is, a union 
of knots. In the case of FS solitons these knots are simple unknots for low Q ) but become 
more complicated as Q grows. For a survey of FS knot solitons, see [14]. 

So the conjecture of [2] results from identifying the FS model as a "submodel" of the 
original TCGL model. In replacing the TCGL model by its FS submodel, one is making a 
two-step truncation. The first truncation, imposing p = p > constant, is reasonable for 
suitable choices of £7, for example 

U = \(pl - - \H 2 ) 2 (5) 

in the limit of large A. Indeed the analogous truncation in ungauged two-condensate systems 
has been found to introduce no drastic qualitative changes [7]. However, the second truncation, 
imposing C = 0, is harder to justify since, even for p constant, C is coupled to (p. Indeed, 
Babaev, Faddeev and Niemi do not directly claim that setting C = is a good approximation 
for the TCGL model [2]. Nonetheless, their claim that TCGL theory should support knot 
solitons assumes that the coupling between and C does not destroy the solitons, and hence, 
that the second truncation C = is, if not a "good approximation", valid at least at a 
qualitative level. This is an assumption, not an established fact, a distinction missed by the 
authors of many of the papers citing [2] , who treat the existence of knot solitons in TCGL 
models likewise as established fact rather than conjecture. 

There have been many attempts to construct knot solitons in TCGL models numerically, 
by choosing ^1^2^ A with tp = [^i,^] having Q > 0, usually in some toroidally symmetric 
ansatz, and minimizing E = J R3 £ using some gradient descent method. In almost all cases 
the initial field configurations have shrunk and fallen through the lattice mesh [10, 11, 19]. 
The only exception we know of is [12], which uses direct gradient flow to evolve the field 
configuration towards a critical point of E ) and claims to have found "definite convergence 
towards torus-shaped configurations" . Given that this finding has not been reproduced by any 
other researchers, it seems likely that this is an artifact of having used too lax a convergence 
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criterion (unfortunately, the precise convergence criterion used is not described in [12], but 
gradient flow is a notoriously slow method, and it is easy to mistake a slowly evolving field 
for a truly static one) . 

The failure of the direct approach to find knot solitons is not surprising given the results 
of [16]. Even if one imposes that p is non- vanishing, so that Q is well-defined, in every degree 
class the infimum of E(i/ji,i/j2,A) is because, for suitably chosen A, any spatially localized 
configuration ^1,^2 is unstable against Derrick scaling. This shows that, for all Q 7^ 0, there 
is no global energy minimizer. Note that this contradicts the common misconception, repeated 
in [12], that the existence of a nontrivial topological charge protects a field configuration from 
collapse. So knot solitons, if they exist in TCGL theory, can only be local minimizers of E, 
and finding them is like looking for a hollow dip in a hill-top, rather than the floor of a valley: 
unless one's initial guess is very close to the dip, one is very unlikely to find it. Indeed this 
point was understood, from a more qualitative, physical viewpoint, by Babaev [1], some years 
after making the original conjecture. For this reason, the failure of previous numerical studies 
to find knot solitons, while discouraging, does not systematically undermine the conjecture of 
[2]. To do so using the direct approach one would have to systematically search the infinite 
dimensional space of all possible initial data. 

In this paper we pursue a different approach. We study the TCGL model in the case 
of a hard-confining potential, (5) with A —> 00, or, equivalently, in the sigma model limit, 
where \ipi\ 2 + |^2| 2 = Po- We may, without loss of generality, take po — 1. This allows us to 
concentrate on the dynamically crucial issue of the coupling between tp and C. Rather than 
minimizing E starting with some plausible initial guess, we consider the one-parameter family 
of models with energies 



parametrized by a G [0, 1], where (•, •) denotes L 2 inner product and || • || denotes L 2 norm. 
Note that at a = 1 this is the TCGL model (in the sigma model limit) while at a = it 
decouples into the FS model and an uncoupled Proca model for C. So the a = model clearly 
supports knot solitons with (7 = 0- they are simply the usual FS knot solitons. The question 
is then, what is the fate of these solitons as the supercurrent coupling is "turned on", that is, 
a is increased to 1? We address this problem numerically by a continuation method. Starting 
at a = 0, we numerically minimize E^ within a given homotopy class using the limited 
memory quasi-Newton algorithm (also called a variable metric algorithm) with BFGS formula 
for Hessian approximations [13]. Having found a minimizer for this value of a ) we increase 
a slightly and minimize E again, starting from the energy minimizer just obtained. In this 
way we construct a curve of energy minimizers parametrized by a ) and the crucial question is 
whether this curve continues all the way to a = 1. 

We will see that, for all a G [0, 1), we have a lower energy bound 



where c > is an absolute constant. This leads one to expect the knot solitons to persist 
whilever a < 1, albeit with nonzero supercurrent C ', and indeed this is what we find. We shall 



E (a) = \\\M\ 2 + 1\\^\\ 2 + f (dc, <p*co) + i||dc|| 2 + Iwcf 



(6) 



E {a) > c VT^\Q\* 



(7) 
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produce convincing numerical evidence, however, that as a —> 1, these knot solitons shrink to 
zero size, because C deforms precisely to the form for which E^ is not stable against Derrick 
scaling. 

Note that our results are not just another negative finding in the same vein as [11], since 
our approach tests not just the conjecture of [2], but also the reasoning underlying it. That is, 
our results conclusively demonstrate that, notwithstanding the formal similarity between the 
TCGL model (in suitable variables) and the FS model, the knot solitons of the latter definitely 
do not persist in the former, even in the most favourable case where a hard confining potential 
enforces \ipi\ 2 + \ip 2 \ 2 = 1- This also elucidates the mechanism by which knot solitons are 
destabilized. It has nothing, in general, to do with p shrinking and aquiring isolated zeroes 
(thus allowing Q to discontinuously drop). Even in the sigma model limit, where Q is a 
rigorously defined topological invariant, the coupling to the supercurrent C alone suffices to 
destabilize the knot solitons. Of course, this does not show for certain that E has no local 
minima with Q ^ 0, but it does invalidate the reasoning used to motivate the conjecture in 
the first place. The present results complement work of Ward [19] (extended by one of us 
[10]) which also embeds TCGL theory into a one-parameter family of models one end of which 
supports solitons. Again, these solitons disappear as the true TCGL theory is approached. It 
also complements the exact results of [16] where it was shown analytically that supercurrent 
coupling destabilizes the Q = 1 Hopf soliton on physical space S 3 . In the absence of evidence 
to the contrary, one should wield Occam's razor and conclude that, in all likelihood, the 
basic TCGL theory does not possess knot solitons. There remains the possibility that more 
elaborate, but physically relevant, versions of TCGL theory, involving direct current-current 
interactions, may have knot solitons [1], but until direct evidence for this is found, we prefer 
to remain sceptical. 



2 Energy bound and Derrick scaling 

We seek local minimizers of the energy functional E^ defined in (6) with Lp : M 3 — >• S 2 and 
C G fi 1 (M 3 ) having the boundary behaviour 

lim ip(x) (foo, lim C(x) (8) 

\x\— >oo \x\— >oo 

and a G [0, 1]. Without loss of generality, we choose cp^ = (0, 0, 1). Configurations with this 
boundary behaviour fall into disjoint homotopy classes labelled by the integer Hopf invariant 



Q = 77T^ / a A da ( 9 ) 

where a is any one- form on M 3 such that da = (p*oj. It is well known that the Faddeev-Skyrme 
energy, which coincides (up to a factor of \ in the usual normalization) with E(°\ip, 0) satisfies 
a topological lower energy bound 

£(°)(^0)>c o |^)|i (10) 
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where c > is an absolute constant [18]. Now, for all a G [0, 1], 

E {a) = ^||cM| 2 + 1(1 - a)||^|| 2 + 1(1 - a)||dC|| 2 + ^\\dC + \<p*u\\ 2 + ^l| 2 
> l||d^|| 2 + l(l-a)||^|| 2 

= ivr^did^f + ii^ii 2 ) (ii) 

where <p(x) = <p(y/l ax). Hence, by the usual Vakulenko-Kapitanski bound (10), 

E( a \<p,C) > coVT^lQ^. (12) 

This leads us to expect that E^ will have a global energy minimizer in each homotopy class 
while ver a G [0, 1). 

It is instructive to subject E^ to the Derrick scaling test for all a G [0, 1] since this gives 
an integral constraint on (<p, C) which we can use as a consistency check on our numerical 
scheme [8]. Assume that (<p, C) is a critical point of E^ a \ Then E^ is stationary with respect 
to all variations of ((/?, C) and so, in particular, with respect to the variation 

y x = <p o ® x , C x = <2)\C (13) 

where A G (0, oo) and Q)\ : M 3 — >• M 3 is the dilation map Qf\{x) = Ax. Note that cpi = and 
C\ — C . A short calculation shows that 

E(X) := ^(^, C x ) = ^||d^|| 2 + ^II^H 2 + ^(<p*u, dC) + ^||dC|| 2 + ^||C|| 2 . (14) 
Hence, since ((/?, C) is a critical point, 

E\l) = -^||d^|| 2 + 1||^|| 2 + dC) + l\\dC\\ 2 - lllCf = 0. (15) 

We shall refer to this as the Derrick constraint. 

Note that, in the case a = 1 (the TCGL model), if C is chosen so that 

dC+^*cj = 0, (16) 
and ((p\, C\) defined as in (13) we have 

^ (1) (^,C A ) = ^||d^|| 2 + ^||C|| 2 ^0 (17) 

as A — V oo. For all ip satisfying (8) there exists C satisfying (16) and (8), which shows that 
EW(<p,C) has infimum in every homotopy class [16]. This argument, which is essentially 
equivalent to the one described in [14] (whose authors attribute it to unpublished work of 
Forgacs and Volkov), explains why the energy bound (12) becomes trivial at a = 1. Of course, 
it does not follow that E^ can have no critical points, because there is no reason why critical 
points of E^ should satisfy (16). 
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For the purposes of numerics, we will study not on IR 3 , but in a large box B = 
[— L, L] x [— L, L] x [— L, L] with Dirichlet boundary conditions (p(x) = y?oo, C(x) = 0, for 
all x G In this case, the Derrick scaling argument must be modified, as follows. First, 

we note that the variation ((p\,C\) is only well defined for A G [l,oo), since for A G (0,1) 
((fx, C\) does not satisfy the boundary conditions. For A > 1 we first notionally extend <p, C 
to the whole of IR 3 by their boundary values, then define ((f\,C\) as in (13). The equation 
(14) for E(X) still holds, but now on [l,oc) rather than (0,oc). We require that E(X) has a 
local minimum at A = 1, but now, since this is at an end-point of the domain, it does not 
follow that E'(l) = 0. Rather, we know only that 

E'(l) = -l\m 2 + i||^|| 2 + ^VdC> + i||dC|| 2 - ^||C|| 2 > 0, (18) 

so the Derrick scaling constraint becomes an integral inequality once we place the theory in a 
finite box. This argument, which could be called "Derrick scaling in a finite box" is of wide 
applicability in classical field theory. Note that it does not really require B C W 1 to be a box: 
any star-shaped domain will do (B C K n is star-shaped if there exists x Q G B such that for all 
x G B the line segment from x Q to x lies in B). Of course, for the purposes of numerics, one 
would hope to choose the box size L to be so large (compared with the soliton core radius) 
that the numerical solution satisfies (15) to a good approximation. We will discuss this issue 
in more detail in section 3. 

3 Numerical results 

We used a simple foward differencing scheme to discretize the energy functional E^ on a cu- 
bic lattice of spacing h with TV 3 points, to give a lattice approximant E^ : (S 2 x IR 3 )^ 3 —> IR. 
Given an initial configuration (<p, C), we used the quasi-Newton BFGS method, as imple- 
mented by the TAO and PETSc parallel numerical libraries [6] and PETSc[3, 4, 5], to find 
a local minimum of E^\ We considered the scheme to have converged to a minimum if the 
sup norm of the gradient of E^ was less than O.Ol/i 3 . We tested the accuracy of this scheme 
by comparing its results at a = 0, for Hopf charges Q = 1, 2, 3, with previous studies of the 
pure FS model [17, 9] and found good agreement (the energy minimizers had the same shape, 
energy and core length to within a few, typically 2, percent). By experimenting in the a = 
case we found that a good balance of accuracy and computational speed was obtained with the 
choices TV = 480 and h = 0.0125, and these values were used for the remaining calculations. 

As described in section 1, having found a minimizer of Eq*\ we incremented a slightly, 
a' = a + 5a, and minimized E^ a ^ starting with the a minimizer as our initial guess. In this 
way, we constructed a curve of minimizers of Eq*\ parametrized by a G [0,1], starting at 
a = and working towards a = 1, in each of the homotopy classes Q — 1, 2, 3. In every case, 
the knot soliton at a = shrinks rapidly as a approaches 1, so that, at a = 1 its core size 
is comparable to the lattice spacing. This extremely coarse configuration can confidently be 
identified as a discretization artifact, rather than a genuine minimizer , as we shall argue 

below. The shrinking process is illustrated in figure 1, which shows a sequence of minimizers 
for increasing a in the case Q — 1. The results for Q = 2, 3 look very similar. A quantitative 
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measure of the shrinking is given in figure 2, which shows the core length of the minimizer 
as a function of a. Recall that the core is by definition the preimage curve of —(foe which, 
in the cases Q — 1,2,3 is a single closed curve. Numerically, we construct this curve using 
interpolation tools of the MayaVi Data Visualizer [15]. This works well until a is very close 
to 1, when the core structure becomes too small for the lattice to properly resolve, and the 
interpolation tool produces a disconnected collection of segments instead of a closed curve. 
From this point, the core length calculation is unreliable, which is why the curves in figure 2 
stop slightly short of a = 1. Clearly, the data up to this point are consistent with the core 
length shrinking rapidly to as a — >• 1. 

Recall we have an energy bound (12) which vanishes at a = 1, and that E^ has infimum 
in every homotopy class, meaning that every class contains fields of arbitrarily low Such 
fields are constructed by Derrick shrinking among the set of fields satisfying (16). Conversely, 
if (16) is not satisfied, the field is protected against shrinking by the presence of a quartic 
term in So if our claim is correct, that the minimizers shrink and vanish as a — >• 1, then 
they should satisfy (16) to a closer and closer approximation as a —> 1. In figure 3 we present 
a graph of ||dC + |(/?*Cc;|| 2 as a function of a which confirms that this is indeed what happens: 
the minimizer loses stability against Derrick scaling as a —> 1, and consequently, its energy 
vanishes rapidly as a — >• 1, as shown in figure 4. Note that these curves are consistent with 
the bound (12). 

Clearly the energy does not vanish exactly at a = 1, but this is a discretization artifact. 
To see this, we define for each a G [0, 1] the quantity 

D{*) = e&I^C) {4 I|CM|2 + l yMl2 + f dC) + ^ l|dC|12 " l 11 ^} (19) 

where (<p, C) is the minimizer of E^ a \ This is the left hand side of the Derrick constraint (15), 
normalized by the total energy of the field. Hence, as shown in section 2, for a minimizer in M 3 , 
one should have D[a) = 0, while for a minimizer in a finite box D[a) > 0. Figure 5 presents a 
plot of D[a) for our numerical minimizers. One sees that D(a) is moderately small for small 
a ) indicating that the boundary of the finite box exerts a moderate but appreciable influence 
on the numerical solution. As a increases D(a) decreases to very close to zero, because the 
minimizers shrink, so the finite box size becomes numerically irrelevant. However, as a gets 
very close to 1, D(a) becomes negative, indicating that from this point on the minimizers 
are strongly affected by discretization effects: if this were a continuum system, any field with 
D(a) < would be unstable against skrinking. In particular, D(l) is large and negative, 
suggesting that the a = 1 "minimizer" is a lattice artifact, and not representative of a genuine 
local minimizer of E^\ 
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(a) a = 



(b) a = 0.7 




(c) a = 0.94 














(d) a: 


= 0.98 



(e) a = 0.99 




(f)a 



Figure 1: Minimizers of for an increasing sequence of values of a, with Hopf charge 
Q — 1. The translucent white surface is the isosurface <p 3 = 0, which can be interpreted as 
bounding the soliton core. The coloured surface is a plot of \C\ on a transverse slice through 
the soliton core. Note that as a —> 1, the core shrinks and the C field accumulates in a singular 
spike. 8 





Figure 3: The squared L 2 norm of AC + \<p*u) of Q = 1 (blue disks), Q = 2 (red squares), and 
Q = 3 (black diamonds) minimizers as a function of a. 
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Figure 4: The energy of the Q = 1 (blue disks), Q = 2 (red squares), and Q = 3 (black 
diamonds) minimizers as a function of a. 
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Figure 5: The Derrick constraint function D(a) of Q = 1 (blue disks), Q = 2 (red squares), 
and Q = 3 (black diamonds) minimizers. This quantity should be non-negative for a genuine 
minimizer in a finite box. 
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